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1. Introduction. 



The dimensionality of space-time plays a key role in all branches of Physics. The 
quantities we calculate, as theoreticians, depend very much on the number of dimensions 
we are considering. Theories in higher and lower dimensions than four have been put 
forth by many researchers and plentiful of good insights have been gained through 
this exercise. Zooming in the arena of quantum field theory, we discover that the 
dimensionality of space-time gained a more sophisticated status, being promoted from 
a mere integer number to that of a complex variable, with the advent and development 
of dimensional regularization by 't Hooft et al and several other pioneers in the field 

In other words, we could say that quantum field theory (QFT), besides other great 
ideas it inspired, physical and mathematical alike, did reveal this amazing possibility: 
the analytic continuation of the space-time dimension D. 

The union between the theory of analytic functions and QFT is very profitable. 
Dimensional regularization (DREG), the technique that bears the concept of analytically 
continued D, is one of its profits. As a step further in this direction Halliday et a/ 0, ^ 
developed the idea of analytically continued D to negative values. Of course, the seminal 
idea of negative values for D is already contained in the work of 't Hooft and others. But 
what is novel in Halliday's insight is the amazing possibility of letting field propagators 
be raised to positive powers, so that the integrand becomes polynomial. The thrust 
behind the idea is that solving a polynomial integral should be — in principle at least 
— easier to perform than rational ones elicited in the usual Feynman integrals. This very 
simple argument, which we call negative dimensional integration method (NDIM), can 
simplify the calculation of Feynman integrals in an astounding way [^, |^, ^, |^, ^, ^ [1^] . 

In the usual DREG 12, |l3l the only quantities that preserve their meaning are 
the Green's functions We will not try to discuss whether they still have (or have 
not) their meaning preserved within the context of NDIM nor speculate what are the 
features, if any, of this "new world" of negative dimensions 0] . What we do is simply to 
allow for it just for calculational purposes. The reader must have this important point 
in mind. 

In our previous works ^ we calculated massive one-loop four point functions 
(former reference) and a massless two-loop three-point vertex (latter reference) with the 
NDIM approach. In the first, NDIM provided not only the well-known hypergeometric 
functions but six other new results in a very straightforward manner [^]; while for the 
two-loop vertex, we considered the particular case where two of its external momenta 
were set on-shell, and NDIM responded with as many as twelve times — surprisingly 
enough, all of them yielding the same correct result — that is, a twelvefold degeneracy. 
This led us to conjecture that when the power series had unit argument and they were 
all summable, then the result would be degenerate. That is, if this conjecture is correct. 
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we need only to carry out one sum — the most convenient one, of course. The conjecture 
remains to be proven or disproved. 

Here we put our NDIM to another "lab-test" by considering two two-loop self- 
energy diagrams which we call by the funny name "flying saucer" diagrams — side view 
(Fig. 1) and front view (Fig. 2), just to make it easier for us to refer to them. The 
outline for this article is as follows: In Section 2 we solve the two-loop Feynman integral 
relative to these two graphs, i.e., the space-time dimension and the exponents of the 
pertinent propagators are left arbitrary. Then, in Section 3, we particularize to suit 
either the "flying saucer, side view" or the "flying saucer, front view" diagram cases. 
And finally. Section 4 is devoted to our concluding remarks. 




Figure 1. Two-loop massless Feynman diagram: the "flying saucer", side view. 



2. Feynman Graphs with Four Massless Propagators. 



The NDIM approach to solve Feynman integrals is beautiful in its simplicity: First, 
we take the propagators of the integral we want to work out, multiply each one of 
them by a specific parameter and then solve the D-dimensional gaussian integral whose 
argument is that very expression. Let us see how it works in practice. Consider the 
gaussian integral. 



J(p^; D) = I d^k d^gexp -aq'-pe --f{p - k)'' ~iu{k - qf 



(1) 



which clearly is relevant to the diagrams we want to work out. It is not difficult to 
integrate it; the result is. 
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Figure 2. Two-loop massless Feynman diagram: tlie "flying saucer" , front view. 

where A = a/3 + Pu + + + auj. Expanding (||) in Taylor series and also expanding 
the multinomial expression in A, we get an eightfold summation, 



oo 



„,=o ni\n2\n:>,\ni\n^\nfi\nj\n^\ 

with the constraint —rii — n2 — n-i — \D = ^4 + ns + ne + + ng coming from the 
multinomial expansion. 

The second step is simpler and faster: expand the exponential (|l|) in Taylor series 
first to get, 

I{p^;D)= E (-1)^^^-^'^-^;^^ (4) 

i,j,l,m=0 v ■ • 

where we define, 

J{z, J, l,m-D) = J d^A; d^g {qJiky [{p - kf]' [{k - qf]"' , (5) 

which is our negative dimensional integral. 

Comparing (3) and (|) we get an expression for the negative-D integral, 

J{z,j,l,m;D) = {-nfip^ycf: , , , ,\ , , (6) 

n~^o ni\n2\n3\n4\n^\neln7lnsl 

where we define the product of gamma functions, 

G = r(i + 2)r(i + j)r(i + /)r(i + m)r(i -a-^D), 
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and since the two expressions must equal, sum indices in the former and exponents of 
propagators in the latter, must satisfy the system, 

ni+n2 + n4 + n6 + ns = i 
ni + ng + n4 + ns = j 
rii + n2 + ris + riQ + rif = I (7) 
n2 + ris + n;s + Tij + ns = m 
rii + n2 + = a 

where a = i+j + l + m + D and the last equation comes from the multinomial expansion. 
Observe that the equations above are linear, but because we have eight unknowns and 
only five equations, in order to solve this system we must choose three of the unknowns 
and solve it in terms of them. There are many ways in which this choice can be done; 
in fact, there are C| = 8!/(5!3!) = 56 possibilities altogether. However, 20 out of the 56 
lead us to trivial solutions, which present no interest at all. The remaining 36 give us 
the results for the Feynman integral when we plug their solutions in equation (P). 

We will solve the non-trivial systems and write down the general results, but before 
doing that, let us see what we can do to lessen our task. Looking at the Feynman 
diagram we can spot symmetry properties that help us in this. Thus, we expect the 
outcoming result to be symmetric under the exchange i ^ m, which in turn will reduce 
by half the number of distinct systems that we need to deal with, since the symmetry 
will account for the remaining half. 

Let us then first consider the solution that leaves ns, uq, ng as free indices in the 
summation; let us call it J'a- It yields. 



{-j - m - |D|r25 - nQ){lD + /[ris + rig) 



where 
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Using one of the properties of the Pochhammer symbol flS 
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one can identify these series as hypergeometric |T6|; in fact, we can rewrite them in a 
convenient manner using another property, 

{a\b + c) = {a\b){a + b\c), (10) 



and sum, for example, the rig series using the well-known formula p5 

T{c)T{c-a-b) 
T{c-a)T{c-b) 



17 ( h \^\ r(c)r(c-a-6) 



yielding, 



(l + l 
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i^D + l\n,)i-l + a\ne) 
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Fs - riQ) 



where 



r(i + z)r(i + /-m) ^^g^ 



T{l-m- lD)T{l + i + l + lD)' 
In a similar manner we can sum the two remaining series, getting as a result 



Ja = i-rcnp 



r(i + z)r(i + /)r(i + m)r(i - a - |d) 



r(l + a)T{l - % - \D)Y{\ -m- \D)V{1 - I - \D) 



T{l+i + 3 +m + \D)T{1 -i-m-D) 
^ V{l + i + m+\D)V{l + l-<y) ' ^ ' 

The last and final step that need to be taken is to bring this result back to our 
real physical world with positive D. Grouping the gamma functions in the numerator 
with the ones in the denominator in convenient Pochhammer symbols and using (P), we 
arrive at 

Jf = vr^(p2)<.(_,|, + ^ + ^£')(-m|i + m + + ]^D\ -2a- ]^D) 

X {-l\a){-i - j - m- ^D\j){i + m + D\-i + l- m- ^D). (15) 

This very simple operation allows us to analytically continue the result back into 
our real physical world, D > 0. Equation (15) is the general result, and we note 
that it is symmetric in i m as it should be, and it is correct [|1^. The reader will 
ask immediately: What is(are) the result (s) that the other solution(s) provide? NDIM 
answers in as brief and surprising a manner as it could be: the same. Indeed just to 
make sure we went through all of them, and verified that it is possible to sum all the 
emerging series and they provide the same result, namely, equation (14) which leads to 
the correct expression, equation (15), i.e., we have a thirty-six-fold degeneracy! 



7 



A word of caution here: Not all the sums can be so easily dealt with. Yet, just 
to convince the reader that it is possible to sum them all, we shall carry out one more 
summation, the hardest one. The degeneracy above mentioned can be classified into 
two sets: 32 solutions are like J a with relatively easy sums to carry out, while 4 of them 
are like the following one which we call jTi,. Consider then the solution with indices 

^1, "-4, ^5, 



ni ,n4, 715=0 



ni!n4!n5!(l - j + a\n4 + n^) 



where 



^ (-J>i+n4 + n5) 
p ^ V{l + i)T{l + l)T{l + m)T{l-a-\D) 

^ r(i-j + a)r(i-m-iD)r(i-2-j-|D)r(i-/-|D)' ^ ' 

Using the same procedure we used for J7a, we can carry out the summation, 
yielding, 

-7 ( ^\D(2.app ^ ik^ + ^I^4)(|I^ + m\m + r24)(-ff |ni) 

i-j\ni + n4) 



X 



where 



(1 - i - j - lD\ni + n^) ' 
T{l-j + a)T{l-l-lD + a) 



T{l + a)T{l-j-l-lD + ay 
However, this time, neither of the remaining sums (in ni or n^) can be written in 
terms of 2-F1, but rather in terms of 3F2 which are not summable for arbitrary values of 
its parameters. Here comes the trick that will do the required job: Put the series in 
terms of a 3F2 function and use the property 0] 

3F2(a,6,c;e, /|1) = Q 3^2(6 - a, / - a, s; s + b,s + c\l), (19) 
where s = e + f — a — b — c and 

_ r(e)r(/)r(.) 
^ r(a)r(s + 6)r(s + c)' 

A good choice is to take 
a = -j + m 
b = \D + l 

c = \D + m + ni (20) 
e = \D + I — a + rii 
f = l-t-3-lD + n, 
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so that the gamma functions in Q simphfy several factors in the series and some of them 
can be grouped by the property (|T0|) giving, 



ni ,714=0 



(lD+j + l-a\ 


\n4){^D + m\ 


\ni){-a\ni) 
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where 



^ (1 - i - ^D\ni){l -i-m-a- Dln^) ^^^^ 
(1 — i + l — m — a — ^D\n4) 



T{^D + I - a)T{l -i-j- ^D)T{1 - i - m - a - D) 



r(-j)r(l -i-a- \D)V{1 -i + l -m- a -\D) 

Now the sum in ni is a 2-^1 function and we can sum it using the Gauss summation 
formula (pJ]). The gamma functions that arise from this summation simplify the 
remaining series in from a 3F2 into a 2F1 function that can be summed again using 
the usual Gauss summation formula. The end result is, 

Jb = {-7r)''{p^rPsP,P5PeP7, (22) 
where Pq came from the rii sum, 

_ r(l - i - cr - lD)T{l -i-m-D) 
^ T{l-i- lD)T{l -i-m-a-Oy 

and P-j came from the last one. 



P7 



r(-j)r(l -i + l-m-a-\D) 



V{\D + l-a)V{l-a + l) ' 

Multiplying all the gamma factors (P3, ■■■jPr) we get, exactly, the expression (14), 
that leads to the correct result (15). 



3. Special Cases. 

The scalar integral we calculated in the previous section has particular cases of interest, 
namely, the "flying saucer" diagrams, side view (Fig.l) and front view (Fig. 2). For the 
side view diagram the exponents of the propagators are all equal to minus one, while 
for the front view diagram the exponents are minus one except for j = —2. 

Let us denote by J'^'^' the general result for the ''flying saucer'^ diagram; when we 
take i = j = / = m = — Iwe have the result for the side view diagram, 

jAC _ 2^D-4 ^^(|-P ' mP " 3)1(2 - 1^)1(4 - D) 

^P> r(3 - mr{D - 2)r(|D - 4) ' ^^^^ 
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while when we take i = —1, j = —2, Z = m = — 1 we have the result for the front view 
diagram, 

jAc _ n(..D-, mD - l)r(5 - D)T{D - 4)r(2 - |D) 



which reproduce well-known results [181 . 



4. Conclusion. 



The analytic continuation of the space-time dimension D into negative values has shown 
us advantages never dreamed of before: we interpret the analytic continuation like 
in usual DREG but solve the Feynman integrals in a much simpler way because the 
integrands we have to deal with are polynomial. The way back road, via another 
analytic continuation is straightforward and the whole procedure is quite simple and 
elegant. There are no cumbersome parametric integrals to solve; on the contrary, the 
only things one needs to know are how to solve gaussian integrals, and systems of linear 
algebraic equations! Furthermore, we have surprising manifold degenerate solutions for 
a single integral. Our previous conjecture on this topic seems to hold, though further 
research is needed to prove it. 
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